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ABSTRACT 

Long-duration gamma-ray bursts (GRBs) are thought to be connected to luminous and energetic 
supcrnovae (SNe), called hypernovae (HNe), resulting from the black- hole (BH) forming collapse of 
massive stars. For recent nearby GRBs 060505 and 060614, however, the expected SNe have not been 
detected. The upper limits to the SN brightness are about 100 times fainter than GRB-associated HNc 
(GRB-HNe), corresponding to the upper limits to the ejected 56 Ni masses of M( 56 Ni) ~ 1O _3 M0. 
SNe with a small amount of 56 Ni ejection are observed as faint Type II SNe. HNe and faint SNe are 
thought to be responsible for the formation of extremely metal-poor (EMP) stars. In this Letter, a 
relativistic jet-induced BH forming explosion of the 40 Mq star is investigated and hydrodynamic and 
nucleosynthetic models are presented. These models can explain both GRB-HNe and GRBs without 
bright SNe in a unified manner. Their connection to EMP stars is also discussed. We suggest that 
GRBs without bright SNe are likely to synthesize M( 56 Ni) ~ 10" 4 to 10~ 3 M@ or ~ 10" 6 M Q . 
Subject headings: Galaxy: halo — gamma rays: bursts — nuclear reactions, nucleosynthesis, abun- 
dances — stars: abundances — stars: Population II — supernovae: general 
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1. INTRODUCTION 

Long-duration gamma-ray bursts (GRBs) at suf- 
ficiently close distances (z < 0.2) have been found 
to be accompanied by luminous and energetic 
Type Ic supernovae (SNe Ic) calle d hypernovae 
(HNe ; GRB 980425/SN 1998bw: iGalama et aT 
19981: GRB 030329/S N 2003dh: iStanek et al 
20031 iHiorth et all l200l GRB 031203/SN 20031w: 
Malesani et all I2004D These GRB-associated HNc 
(GRB-HNe) are suggested to be the outcome of very 



energetic bl ack hole (BH) forming 
stars (e.g., Ilwamoto et al.l I1998T 



3 Ni) 
4 - 



explosions of massive 
The ejected 56 Ni 
~ 0.3 - 0.7M Q (e.g., 
10 times larger than 



O.O7M ; iNomoto et all 

nearby long-du ration 
Fvnbo et all [2006) and 



mass is estimated t o be M( 56 
iMazzali et"aT]|2006ft which is 
typica l SNe Ic [M( 56 Ni) ~ 
I2006bj . 

For recently discovered 
GRB 060505 (z = 0.089, 

GRB 060614 (z = 0.125. IGal- Yam et al.l 120061: 
Fvnbo et all 120061 : iDella Valle et al J 120061 : iGehrels et all 
2006|) . no SN was detected. 7 Upper limits to brightness 
of the possible SNe are about 100 times fainter than 
SN 1998bw [i.e., M( 56 Ni) < 10- 3 M©]. 
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A small amount of 56 Ni ejection has been indi- 
cated in the fain tness of several Type II SNe (SNe II, 
e.g., SN 1994W. iSollerman et al.l fl998t and SN 1997D, 
ITuratto et al.l ll998[ L The estimated explosion ener- 
gies, E, of these faint SNe II are very small (E < 

10 51 ergs, ITuratto et al.1 I1998D . These properties are 
well-reproduced by the spherical explosion models that 
undergo significant fallback if E is sufficient ly small 
(jWooslev fc Weave j [1991 Ilwamoto et al.l I2005D . Thus 
these faint SN explosions with low E seem to be su- 
perficia lly irreconcilable to the formation of energetic 
GRBs (IGal- Yam et al.l l2006f). However. Nomoto et al. 
(2006a: 8 see also lNomoto et al.ll2004D has predicted the 
existence of "dark hypernovae" (i.e., long GRBs with no 
SNe) based on the argument in the next paragraph. In 
this Letter, we present actual hydrodynamical models in 
which a high-energy narrow jet produces a GRB and a 
faint/dark SN with little 56 Ni ejection. 

An indication that some faint SNe produce high- 
energy jets is seen in the abundance patterns of 
the extremely metal-poor (EMP) stars with [Fe/H] 
< —3. 5. 9 It has been suggested that the abun- 
dance patterns of these EMP stars show the nucle- 
osynthesis yields of a single core-collapse SN (e.g., 
iBeers fc Christliebl .2005). In particular, the C-enhanced 
type of th e EMP stars have been well explained by the 
faint SNe (lUmeda fc Nomotoll2005L Ilwamoto et~aTll2005l : 
INomoto et al.l l2006bt iTominaga et al.l l2007h. except for 
their large Co/Fe and Zn/Fe rat ios (e.g.. iDepagne et al.l 
I2002t iBessell fc Christliebl 120051) . The enhancement of 
Co and Zn in low metallicity stars requires explosive nu- 
cleosynthesis under high entropy. In a "spherical" model, 



8 The ppt file is available from "program" at http://www.me- 
rate.mi.astro.it /docM/OAB /Research/SWIFT /sanservolo2006 / 

9 Here [A/B] = log 10 (Ar A /jV B ) - ^g w (N A /N B ) Q , where the 
subscript refers to the solar value and 7V A and Nb are the abun- 
dances of elements A and B, respectively. 
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a high-entropy explosion corresponds to a high-energy 
explosion that inevitably synthesizes a large amount 
of 56 Ni. One possible solution to this incompatibil- 
ity is that some faint SNe are associated with a nar- 
row jet within which a h igh-entropy region is confined 
(jUmeda fe Nomotol 120051 ). If this would be a realistic 
model for the EMP stars, some faint SNe would accom- 
pany sufficiently energetic jets to produce GRBs. 

In this Letter, we present hydrodynamical and nuclc- 
osynthetic models of the 40 M star explosions with rela- 
tivistic jets. 10 We show that these models can explain the 
existence of GRB-HNe and GRBs without bright SNe in 
a unified manner. We suggest that both GRB-HNe and 
GRBs without bright SNe are BH-forming SNe involving 
relativistic jets and that these explosions are responsible 
for the formation of the EMP stars. 

2. MODELS 

We investigate the j et-induced exp losions (e.g., 
iMaeda fc Nomoto I l2003t INagataki et~all 120061) of the 
40M^ Population II I stars (jUmeda fc Nomotol 120051 : 
iTominaga et al.ll2007f ) using a multi-dimen sional special 
relati vistic Eulerian hydrodynamic code (jUmeda et al.l 
120051: N. Tominaga, 2007, PhD thesis, in preparation). 
We assume that the explosions of the 4OM Popula- 
tion III stars involve relativistic jets in analogy with the 
4OM Pop I GRB-HNe, since the Fe core masses of Pop- 
ulation III and Population I stars are similar. 

We inject the jets at a radius R[ n ~ 900 km, corre- 
sponding to an enclosed mass of M ~ 1.4M , and follow 
the jet propagation. Since the explosion mechanism of 
GRB-HNe is still under debate, the jets are treated para- 
metrically with the following five parameters: energy de- 
position rate (Pdep), total deposited energy (Pdcp); ini- 
tial half angle of the jets (#j t), initial Lorentz factor 
(Tjot), and the ratio of thermal to total deposited ener- 
gies (/th). 

In this Letter, we investigate the dependence of nu- 
cleosynthesis outcome on E^ ep for a range of Pdcp,5i = 
^dcp/lO^ergss- 1 = 0.3 - 1500. The diversity of 

.Edep is consistent with the wide range of the observed 
isotrop ic equivalent 7-ra y energies and timescales of 
GRBs ()Amati et al.ll2006l and refer ences therein). Vari- 
ations of activities of the central engines, possibly cor- 
responding to different rotational velocities or magnetic 
fields, may well produce the variation of Pdep- We expe- 
diently fix the other parameters as Pdep = 1-5 X 10 52 ergs, 
0jet = 15°, r jct = 100, and f th = 10~ 3 in all models. 

The thermodynamic history is traced with marker par- 
ticles representing indivi dual Lagrangian elements (e.g., 
IMaeda fc Nomotol 120031 ) . The nucleosynthesis calcula- 
tion is performed as pos t-processing and includes 280 
species up to 79 Br (see lUmeda fc Nomotol I2005L Ta- 
ble 1 ). Because rec ent studies (e.g.. iRampp fc Jankal 
120001 iFrohlich et al.lf2006h suggest that the electron frac- 
tion (Y c ) is affected by neutrino processes, we assume 
Y e = 0.5001 in the Fe cor e to maximize the Co/Fe ratio 
(lUmeda fc Nomotol l2005h . 

The thermodynamic history of the jet is also traced 

10 Details of the numerical method, input physics, and the per- 
formance of the code will be described in N. Tominaga (2007, PhD 
thesis, in preparation). 



by the marker particles. In computing the jet composi- 
tion, we assume that the jet expands adiabatically from 
the Schwartzshild radius to the inner boundary of the 
computational domain, and that the jet initially has the 
composition of the accreted stellar materials. The details 
will be presented in future work. 

3. RESULTS 

In order to inject the jet, the ram pressure of the 
jet (Pjot) should overcome that of the infalling material 

(-Pfaii)- Pjet is determined by R in , E dcp , Ojet, Pjet, and / th , 
thus being constant in time in the present models. On the 
other hand, Pf a n decreases with time, since the density of 
the outer materials decreases following the gravitational 
collapse (e.g.. iFrver fc M cszaros 2003). For lower Pdep, 
Pjet is lower, so that the jet injection {Pj C t > Pali) is re- 
alized at a later time when the central remnant becomes 
more massive due to more infall. 

After the jet injection is initiated, the shock fronts be- 
tween the jets and the infalling material proceed outward 
in the stellar mantle. Figure[T]is a snapshot of the model 
with Pdcp,5i = 15 at 1 s after the start of jet injection. 
When the jet injection ends, the jets have been decel- 
erated by collisions with the dense stellar mantle and 
the shock has become more spherical. The inner mate- 
rial is ejected along the jet-axis but not along the equa- 
torial plane. On the other hand, the outer material is 
ejected even along the equatorial plane because the in- 
fall along the equatorial plane is ter minated as the shock 
reaches the equatorial pla ne (e.g., IMaeda fc Nomotol 
120031: INagataki et al.ll200ffl . 

3.1. 56 Ni Mass 

The top panel of Figure [2] shows the dependence of 
the ejected M( 56 Ni) on the energy deposition rate Pdep- 
For lower Pdop, smaller M( 56 Ni) is synthesized in ex- 
plosive nucleosynthesis becaus e of lower post-shock den- 
sities and temp eratures (e.g.. IMaeda fc Nomoto I [20031 : 
INagataki et all 120061 : Maeda & Tominaga 2007, here- 
after IMT07D . While the materials in the C+O layer fall 
through the inner boundary, Pf a n decreases only mod- 
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dyn cm 2 ) during this phase because of the relatively 
shallow pre-SN density gradien t in the C+O layer (e.g., 
IFrver fc Meszarosll2003l : iMTOl . 

If Pdcp,5i > 3, Pjet > Pf^ so that the jet injection is 
initiated below the bottom of the C+O layer, leading to 
the synthesis of M( 56 Ni) > lQ- 3 M e . If P dop ,5i < 3, on 

the other hand, Pj ot < P^ ° so that the jet injection is 
delayed until Pf a n decreases below Pj et and initiated near 
the surface of the C+O core; then the ejected 56 Ni is as 
small as M( 56 Ni) < 1O- 3 M . 

56 Ni contained in the relativistic jets is only M( 56 Ni) ~ 
10~ 6 to 1O _4 M because the total mass of the jets is 
Afjot ~ 1O~ 4 M in our model with Pj et = 100 and 
Pdcp = 1.5 x 10 52 ergs. The 56 Ni production in the jets is 
predominant only for Pdep,5i < 1-5 in the present model 
(Fig. [2]). The models cannot synthesize 56 Ni explosively 
and eject very little M( 56 Ni) (~ 1O- 6 M ). 

3.1.1. GRB-HNe 
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For high energy deposition rates (-Edep, 51 > 60), the 
explosions synthesize large M( 56 Ni) (> 0.1M Q ) being 
consistent with GRB-HNe. The remnant mass was 
M ~ 1.4M Q when the jet injection was started, but it 
grows as material is accreted from the equatorial plane. 
The final BH masses range from Mbh = 10.8Mq for 

Sdo P ,5i = 60 to Mbh = 5.5M Q for £ d ep, 5 i = 1500, 
which are consisten t with the masses of stellar-mass BHs 
(jBailvn et al.lll998l ). The model with £ do p, 5 i = 300 syn- 
thesizes M( 56 Ni) ~ 0.4M Q and results in M B h = 6.4M . 

Since the jet injection with large Ed cp is in short 
timescale and terminated before the jets reach the sur- 
face of the C+O core, the asphericity of the ejecta inside 
the C+O core for Edep, 51 > 60 is probably consistent 
with a relatively oval explosion as indicated by the l ight 
curve and spectra of GRB-HNe (jMaeda et al.N2006aH bl). 

Neutrino annihilation has been estimated to provide 
on ly a small ene rgy deposition rate (e.g., Edep, 51 ~ 
1; IWooslevi 1 19931 ). leading to the synthesis of smaller 
M( 56 Ni) than is required for GRB-HNe. In order to 
synthesize a sufficient amount of 56 Ni, large -Edep should 
be produced via som e mechanism, e.g., magneto- rotation 
(M izuno et afll2004t another po ssibility , 56 Ni production 
in the disk wind, is discussed in IMT07I ). 

3.1.2. GRBs 060505 and 060614 

For low energy deposition rates ( -Edep, 51 < 3), the 
ejected 56 Ni masses [M( 56 Ni) < 10 _3 M Q ] are smaller 
than the upper limits for GRBs 060505 and 060614. 
The final BH masses range from Mbh = 18.2M Q for 
£de P ,5i = 0.3 to Mbh = 15.1M for £ dcp ,5i = 3. 
While the material ejected along the jet direction in- 
volves those from the C+O core, the material along the 
equatorial plane is ejected only from the outer part of 
the H envelope. Thus Mbh exceeds the C+O core mass 
M c +o = 13.9M Q for the 40M Q star. 

If the star lost the H and He envelopes before its 
core collapsed, and if the explosion is viewed from the 
jet direction, we would observe GRB without SN re- 
brightening. This may be the situation for GRBs 060505 
and 060614. 

3.1.3. GRBs with Faint or Subluminous SNe 
For intermediate energy deposition rates (3 < 
-Edep, 51 < 60), the explosions eject 10~ 3 M Q < 
M( 56 Ni) < 0.1M Q and the final BH masses are 
10.8M© < Mbh < 15.1M Q . The resulting SN is 
faint [M( 56 Ni) < 0.01M©] or sub-luminous [0.01M Q < 
M( 56 Ni) < O.1M ]. 

Nearby GRBs with faint or sub-luminous SNe have not 
been observed. Possible reasons may be that (1) they do 
not occur intrinsically, i.e., the energy deposition rate is 
bimodally distributed, or that (2) the number of observed 
nearby GRBs is still too small. For case 1, the GRB pro- 
genitors may be divided into two groups, e.g., with rapid 
or slow rotation and/or with strong or weak magnetic 
fields. For case 2, future observations will detect GRBs 
associated with a faint or sub-luminous SN. 

3.2. Abundance Ratio: C/Fe 



The bottom panel of Figure [2] shows the dependence of 
the abundance ratio [C/Fe] on -Edep- Lower Edep yields 
larger Mbh and thus larger [C/Fe], because the infall 
decreases the amount of inner core m aterial (Fe) relative 
to th at of outer material (C) (see also lMaeda fe Nomoto I 
|200"1 . As in the case of M( 5 %i) (§3.1), [C/Fe] changes 
dramatically at Edep, 51 ~ 3. 

The abundance patterns of the EMP stars are good 
indicators of nucleosynthesis in a single SN because the 
Galaxy was effe ctively unmixed at [Fe/H] < —3 (e.g., 
lTumlinsonll2006h . They are classified into three groups 
according to [C/Fe]: (1) [C/ Fe] ~ 0, normal EMP 
stars (-4 < [Fe/H] < -3, e.g.. iCavrel etaHl200"l : (2) 
[C/Fe] > +1, Carbon-enhanced EMP (CEMP) stars 
(-4 < [Fe/H] < -3, e.g., CS 22949-37: iDeoagne et all 
120021 ): (3) [C/Fe] ~ +4, hyper m etal-poor (HMP) stars 
([Fe/H] < -5, e.g., H E 0107-5240:lCh ristlieb et al.ll200l 
Bessellfe Christliebl [20051 HE 1327-2326. iFrebel et all 
2005). 

Figure 3 shows that the general abundance patterns 
of the normal EMP stars, the CEMP star CS 22949- 
37, and the HMP stars HE 0107-5240 and HE 1327- 
2326 are reproduced by models with Edep, 51 = 120, 3.0, 
1.5, and 0.5, respectively (see Table 1 for model pa- 
rameters). The model for the normal EMP stars ejects 
M( 56 Ni) ~ 0.2M Q , i.e. a factor of 2 less than SN 1998bw. 
On the other hand, the models for the CEMP and the 
HMP stars eject M( 56 Ni) - 8 x 10~ 4 and 4 x 10~ 6 M Q , 
respectively, which are always smaller than the upper 
limits for GRBs 060505 and 060614. The lack of the 
metal-poor stars at — 5 < [Fe/H] < —4 is explained by 
the narrow range of Edep- The N/C ratio in the models 
for CS 22949-37 and HE 1327-2326 is enhanced by par- 
tial mixin g between the He and H layers during pre-SN 
evolution (jlwamoto et al.ll2005l ) . 

4. DISCUSSION AND CONCLUSION 

We have computed hydrodynamics and nucleosynthe- 
sis for the explosions induced by relativistic jets. We 
have shown that (1) the explosions with large Edep are 
observed as GRB-HNe and their yields explain the abun- 
dances of normal EMP stars, and (2) the explosions with 
small E dcp are observed as GRBs without bright SNe and 
are responsible for the formation of the CEMP and the 
HMP stars. We thus propose that GRB-HNe and GRBs 
without bright SNe belong to a continuous series of BH- 
forming SNe with the relativistic jets of different Edcp- 

Presently, the number fraction of GRBs without bright 
SNe relative to the know n nearby GRBs is ~ 40% (e.g., 
IWooslev fe Blooml l2006h . On the other hand, among 
the EMP stars with [Fe/H] < -3.5, the fractions of the 
CEMP and the HMP stars relative to the EMP stars 
are ~ 25% and ~ 15%, respectively dBeers fe Christliebl 
2005). Although the numbers of observed GRBs and 
the EMP stars with [Fe/H] < -3.5 are still too small to 
discuss statistics, the fraction of GRBs without bright 
SNe is consistent with the sum of the fractions of the 
CEMP and the HMP stars. Thus GRBs 060505 and 
060614 are likely related to the CEMP stars ejecting 
M( 56 Ni) ~ 10~ 4 to 1O" 3 M or HMP stars ejecting 
M( 56 Ni) - 10- 6 M Q . 

A short GRB , probabl y the result of the m erger of two 
compact objects (e.g.. iGehrels et al"1l2005l ). synthesizes 
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TABLE 1 

Models compared with metal-poor stars. 



Stars 


^dcp 


M( 56 Ni) 




[C/Fc] 




[10 51 ergs s- 1 ] 


[M e ] 


[M @ ] 




EMP 


120 


2.09 x 10" 1 


9.1 


0.02 


CS 22949-37 


3.0 


7.64 x 10~ 4 


15.1 


1.20 


HE 1327-2326 


1.5 


3.90 x 10~ e 


16.9 


3.21 


HE 0107-5240 


0.5 


2.80 x 10~ 6 


17.1 


3.91 



virtually no 56 Ni because the ejecta must be too neutron- 
rich. In contrast, our model suggests that GRBs 060505 
and 060614 produced M( 56 Ni) ~ 10~ 4 - 10" 3 Af Q or ~ 
10~ 6 A/ Q . If such a GRB without a bright SN occurs 



in a very faint and nearby galaxy, our model predicts 
that some re-brightening due to the 56 Ni decay can be 
observed. 
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Fig. 1. — Density structure of the 40 Mq Population III star explosion model of £?dep,5l = 15 at 1 s after the start of the jet injection. The 
jets penetrate the stellar mantle (red arrows) and material falls on the plane perpendicular to the jets (black arrows). The dots represent 
ejected Lagrangian elements dominated by Fe ( 56 Ni, red) and by O (blue). 
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Fig. 2.— Top: Ejected 56 Ni mass (red: explosive nucleosynthesis products, blue: the jet contribution) as a function of the energy 
deposition rate. The background color shows the corresponding SNe (red: GRB-HNe, yellow: sub-luminous SNe, blue: faint SNe, green: 
GRBs 060505 and 060614). Vertical lines divide the resulting SNe according to their brightness. Bottom: Dependence of abundance ratio 
[C/Fe] on the energy deposition rate. The background color shows the corresponding metal-poor stars (yellow: EMP, red: CEMP, blue: 
HMP stars). 
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FlG. 3. — Comparison of the abundance patterns of metal-poor stars and of models. Top: Normal EMP (red dots) and CEMP (blue 
triangles) stars and models with -Edep,5i — 

120 (solid line) and = 3.0 (dashed line). Bottom: HMP stars: HE 1327-2326, (red dots), and 
HE 0107—5240, (blue triangles) and models with £?dcp,5l = 1.5 (solid line) and = 0.5 (dashed line). 



